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The new numerical approach for consideration of quantum dynamics and calculations of the 
average values of quantum operators and time correlation functions in the Wigner representation 
of quantum statistical mechanics has been developed. The time correlation functions have been 
presented in the form of the integral of the Weyl's symbol of considered operators and the Fourier 
transform of the product of matrix elements of the dynamic propagators. For the last function the 
integral Wigner- Liouville's type equation has been derived. The initial condition for this equation 
has been obtained in the form of the Fourier transform of the Wiener path integral representation 
of the matrix elements of the propagators at initial time. The numerical procedure for solving this 
equation combining both molecular dynamics and Monte Carlo methods has been developed. An 
application of the developed approach to the micro canonical ensemble has been also considered in 
the second part of this paper. 

For electrons in disordered systems of scatterers the numerical results have been obtained for 
series of the average values of the quantum operators including position and momentum dispersions, 
average energy, energy distribution function as well as for the frequency dependencies of tensor of 
electron conductivity and permittivity according to quantum Kubo formula. Zero or very small 
value of static conductivity have been considered as the manifestation of Anderson localization of 
electrons in ID case. Independent evidence of Anderson localization comes from the behaviour 
of the calculated time dependence of position dispersion. Nevertheless for localized electrons the 
energy distribution function obtained has the long exponentially decaying tail, which is the reason 
of exponentially rare appearance of large values of quantum particle virtual energy that strongly 
affects the behaviour of the position dispersion. 
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I. INTRODUCTION 



It is known that crystalline materials have been studied intensively by physicists for a long time since the beginning 
of quantum physics. Practically the same is valid for investigations of the disorder state of the matter. Disorder 
exist in varying degree, ranging from the a few impurities or interstitials in an otherwise perfect crystalline host to 
the strongly disordered limit of alloys or glassy structures. In the past few years there has been growing realization 
that for understanding the disordered materials new concepts must be introduced which treat the disorder from the 
very beginning E[. The new understanding is based on advances in two different areas . The first is the problem 
of Anderson localization, which deals with the nature of the wave function of single electron in the presence of the 
random potential. A scaling description of the Anderson localization problem has greatly deepened understanding 
of the problem. The second aspect of the problem is the interaction among electrons in the presence of a random 
potential. The fact that the electrons are diffusive instead of freely propagating leads to a profound modification of 
the traditional view based on the Fermi-liquid theory of metals. 

In this paper we treat the first of the mentioned above problem. In 1958, Anderson pointed out that if disorder is 
strong enough, the electron wave function may become localized , in that the its envelope decays exponentially from 
the initial electron position in space. The characteristic space scale of this decay is called as localization length. The 
energy and density dependence of localization length results in the metal- insulator transition and the concepts of the 
mobility energetic edge. Despite of significant advances the mentioned and a lot of other related problem are unsolved 
now due to thire complexity. For investigation of the problem of Anderson localization we have applied the developed 
in , H numerical approach allowing rigorous calculations of the quantum time correlation functions and average 
values of quantum operators in the Wigner formulation of quantum statistical mechanics. 

Within Wigner formulation of quantum mechanics the time correlation function can be presented in the form of the 
integral of Weil's symbols of the operators F and A and so called spectral density, which is the Fourier transform of 
the product of the dynamic propagator matrix elements arising due to Heisenberg representation of the operator A(t). 
For spectral density the integral equation of the Wigner- Liouvillc' type has been derived. The developed approach 
for solving this equation combines both molecular dynamics and Monte Carlo methods and allows to combine the 
existing molecular dynamics and Monte Carlo codes after small modernization. Efficiency of the developed method 
is high enough to use for calculations not very fast computers. The explicit expression of the initial condition for 
this integral equation is convenient to obtain as Fourier transform of the propagator matrix elements presented in the 
form of Wiener path integrals. 

The numerical results have been obtained for series of the average value of the quantum operators as well as for the 
frequency dependencies of tensor of electron conductivity and permittivity according to quantum Kubo formulas for 
electrons in disordered systems of scatterers. Zero or very small value of static conductivity have been obtained as the 
manifestation of Anderson localization of electrons in ID case. Independent evidence of Anderson localization comes 
from the behaviour of the calculated time dependence of position dispersion. Nevertheless for localized electrons the 
energy distribution function obtained has the long exponentially decaying tail, which is the reason of exponentially 
rare appearance of large values of quantum particle virtual energy that strongly affects the behaviour of the position 
dispersion. 



II. WIGNER REPRESENTATION OF QUANTUM STATISTICAL MECHANICS 



The time correlation functions Cf a(£) = ( F(0)A(t)^ for different dynamic properties are among the most inter- 
esting quantities in studying dynamics of electrons in disordered systems of scatters, transport coefficients, chemical 
reaction rates, consideration of equilibrium and transient spectroscopy and so on. Our starting point is the general 
operator expression for canonical ensemble averaged time correlation functions B : 



C FA (t) = Z^Tr (Fcxp (iHt*Jhj icxp (-iHt c /tiJ 



Here H is the Hamiltonian of the system, H = K + U, K is the kinetic energy operator, U is the potential energy 
operator, t c = t — ihf3/2, (3 — l/fc^T, F and A are quantum operators of the considered dynamic quantities, 

Z = Tr ^exp (—/3HJ ^ is the partition function. Wigner representation of the time correlation function can be 



written as: 



Cfa (*) = 1 1 dp 1 dqidp 2 dq 2 F (pi, <7i) A (p 2 ,<?2) 

xW (pi ! qi;p2,q2;t;ih/3) 
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where the spectral density W {pi,qi;p 2 ,q2',t; iU(3) is defined by: 



Wipx^p^tiihp) = Z^JJdCid^ exp (i^j exp {i^f) 



x + T 



exp 



g 2 -4f)(g 2 + 4f 



exp I —iHt c /h 
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and F (pi, gi) and A (^2,^2) are Weyl's symbols of operators H F and A : 



A (p, g) 



d£exp 



where (g' 
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q"j are the matrix elements and v is the space dimension. So the problem of numerical calculation of the 

canonically averaged time correlation function can be reduced to the consideration of the spectral density evolution 
describing as can been proved by the following integral equation : 



W(pi,qi;p2,q2;t;ihf3) = W(po,qo;po,q ;ihf3)+ 
J Q dr J dsdi]W(p T - s,q T ;p T - r\, q T ; r; ih(3)~j(s, q T ; r), gy) 



(1) 



where j(s, g r ; 77, q T ) = h{w (s, q T ) 5{rf)~w (77, q T ) <5(s)}, 6(s) is Dirac delta function, w (s, g) is defined by the expression 



w(s,q) 



d q >U(q-q') S m(^)+F(q) 



dS (s) 
ds 



F (q) is the classical force, {q T (r;pi, qi,t), p T (r;pi, qi,t)} and {q T (r;p2, q2,t), p T (r\pi-, g 2 , t)} are the pare of dynamic 
pg-trajectories for 'negative and positive time direction' respectively and initial condition at r = t: 



dp/dr = F(g T (r))/2; q t (t;pi,qi,t) = qi 
dq/dr = p T (r)/2ro; Pt(t;Pi,qi,t) =pi 



(2) 



dp/dr = -F(q T (r))/2; q t (t;p 2 ,q2,t) = q 2 
dq/dr = -p T (r) /2m; p t (t;p 2 , ?2,*) = P2 

The spectral density initial condition W(pi t qi;p2,q2', ihjS) = W (pi,gi;p2,g2,; 0; i%(3) can be written in the form of the 
finite difference approximation of the path integral : 

W (p 1 ,q 1 ;p2,q2;ih/3) » 
« / / dqi...dq M J J dq[ „dq' M V (pi, gi;p 2 , 92; 9l, -,9m! 9i) ■■■^'m^^-P) 
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exp I — eF 



*(pii 915^2, ga; 91, -,9m; gl, -,<i'u' ih 0) 

exp ( — eF 

.(-*) 



exp ^— eF^j 

...exp(-et/(9Af)) (gM 



gi^)exp(-eC/ (gi)) (gi 



exp 1 



gi)exp(-e?7(gi))(g; 



-exp(-eC/ (<4f)) 
^ (p;g',g") = (2A 2 ) 1,72 exp f - 



g 2 )exp(-eZ7(g* 2 )) x 



92) <p(P2;qM,q[) x 



exp I — eF 



exp ( — eF 



g 2 )exp(-e£7 (g 2 )) x 



9i)y(Pi;9M,9i) 



(pA/R+»7r(g'-g")/A| pA/M-wr^ -g")/x) 
2tt 



(3) 



where e = ,M»1 and A 2 = 2irh 2 l3/2mM . 

Let us rewrite the integral equation (0) and the iteration form of its solution in symbolic form: W l = W* + K^W 7 
and 



W l — W l + K l Ti W T1 + K* 2 K%W Tl + K* a K%K™W Tl 



(4) 



Here W l and W Tl is the quantum initial density evolving classically in intervals [0,i] and [0,ri], while Ft* +1 are 
operators, which describe propagation between times 7j and Tj+i. The time correlation functions are the linear 
functionals of the spectral density: 
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C F A (*) = (2^)^ / / dpidq 1 dp 2 dq 2 F {p u qi) A(jp 2 M2) X 

xW (pi,qi;p2,q2;t;ihP) = (5) 
(<f>\W*) + (flK^W 7 *) + {4\K^K^W n ) + (flK^KZKgW 7 *) + ... 

where brackets (|) for functions <f> (pi,gi;p2,?2) = F (pi, gi) A (pa.Qte) and W(po> <7o;Po, go; or K*.K%_ 1 ...K^W T1 
mean the integration over the phase spaces {p± t q\\ P2,Q2 }• 

Note that average values of quantum operators i* 1 (t) can be formally presented in the form analogous to (|5|) : 

F(t) = Z~ x Tr (exp fiHt/h \ Fexp f-iHt/hj p(0) \ = 

- ]_,2v I I dpidqidp 2 dq 2 -{F(pi,qi)+F(p2,q2)}x 
(27m) J J ^ 

xW {p 1 ,qi\p2,q 2 \t]ih(3) 

where p (0) is the initial density matrix and Z = Tr (p (0)). To obtain the explicit expression of the terms of series (^) 
and to analyze its mathematical structure let us introduce the pieces of dynamic trajectories defined by (^|) according 
to the following recurrent relations: 



p 

q 



-I =3+1 =3+1 \ 

-1 =P{T-j-l,Pj - S ji1j > r j) 

-/ =3+1 =3+1 ^ 

~i -j+i =3+1 ^ W 

-1=^-1;^ 



-/ -7+1 \ 

for initial conditions at times Tj satisfying the following inequalities {t = Tj > ... > T2 > T\ > tq = 0, j = 0, 1, ...}. If 
Tj =fwe assume that Sj = 0, p^ +1 = pi t gj +1 = q\ and 77^ = 0, = £>2j = <?2 and superscript in the left part 

of the expression is omitted pj_i = Pj-i, Qj-i — Qj-i and pj_x = Pj-i, Qj-i — Qj-i- So for example the third term 
(j = 3) can be written as: 



t PT3 

dr 2 \ 
11 Jo 



KLK^W 1 = dr 2 dn / ds 2 d7 l2 / d Sl d m x 



xj(s 2 , q2;V2,q2h(si, ql;m, qiW(pl,ql;p~l, ql;ih(3) 



III. QUANTUM DYNAMICS 

The possibility to convert series like (|^) into the form convenient for probabilistic interpretation allow us to develop 
the Monte Carlo method for its calculation 0] , || . According to the general theory of the Monte Carlo methods for 
solving linear integral equations one can simultaneously calculate all terms in the iteration series like (Q) . So using 
the basic ideas of M we have developed the Monte Carlo scheme , which provides domain sampling of the terms giving 
the main contribution to the series (||) . This sampling reduces also numerical expenses for calculations of integrals of 
each terms . Ensemble averaging on the configuration of interacting quantum particles and classical heavy scatterers 
has been performed according to the probability distribution proportional to (^), while the dynamic evolution has 
been realized according to equations (^|) . 

IV. SYSTEM OF UNITS. 

Numerical calculations are more convenient to perform for dimensionless equations. Let the quantum system 
considered have the following characteristic energy Eq and time To scales. The dimensionless combination EoTq/K 
may be considered as a measure of quantum behaviour of a system. However in our quantum dynamics studies it is 
more convenient to fix the maximum value of the considered time t' and to take it as a unit of dynamic time t of the 
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system. So dimensionless time r = t/t' will always vary from up to 1. As unit of a length we take the reciprocal 
wavenumber fc _1 , determined by the ratio k 2 — 2mEo/h 2 , where m is characteristic mass of quantum particle. So, 
for example, for one electron in the external potential field VqU (q) the operator exponent of the time propagator can 
be rewritten in the form: 

Ht/h = I -—A + V U (q) i t/h = 

■ — A + — —U(q)\T = I — A + £ U(q)\T 

2M 2E M W J \ 2M M W J 

Here A is Laplace operator, Vq is characteristic constant of interaction of potential field, M = h/2Eot', £o = 
Vq/2EoM = Vot' /h, t = t/t'. The similar expressions and system of units were made use of for one electron in a field 
of chaotic classical scatterers. In this system of units all the above mentioned formulas take more simple form as h 
can be substituted for by the unity, physical mass m by formal parameter M and constant of interaction is redefined 
by relation £ = V^t' /h. 

A computation of the tensor of the electrical conductivity has been performed as the first example of application of 
the developed approach. According to quantum Kubo formulas in one electron approximation the tensor of electrical 
conductivity is 



(oj) = n j exp (iut — et) J (^J 7 J a (t + ih\)^)d\dt 



where e — > 0, n is the electron density, J a — eq a (t) = ep a /m is the electrical current operator , q a is the a component 
of electron velocity operator. Wigner representation of this tensor may be written in the form: 

<7q, 7 (cj) = n J °° exp [iiot — et) (j) aj (t, A) dXdt = 
= nk~ 2 e 2 o al (oj) /2Eot' 

lp a -y (t, A) = (^J-yJa (t + ih\)^ = 

= (2 X,f" I dpidqidp 2 dq 2 J 7 (pi,qi) J a (p2,q2)W{pi,qi;p 2 ,q2;t;ihP;ih\) 

W(p 1 ,q 1 ; P 2, q2;t; ifi(3- ifi\) = J d^d& exp (i^ exp (i^) 
9i + I | exp (i*f) | 92 - %) U + % |exp {-i^f) \ Ql - | 



where ri = ti + ih\„ t 2 = t 2 — ih((3 — A). Here a ai (to) is the dimensionless tensor of electrical conductivity. In 
the model of independent electrons in the medium of chaotic classical scatterers (one particle approximation) the 
Hamiltonian of the system can be written in the form: 

H = H es + Hg S 

N 



H es = -A/2M + J2to s U es (\q - QA/*') 



N N 

Hf B = Y,p 2 ™/2M™s + Y,Zo s Uss (\Qi - Qj\/o") 

3 = 1 i^j 

where m s is the mass of scatterer (m/m s -C 1), Hf s is the classical Hamiltonian of scatterers, £g s = V^t' /Ti and 
Co" — VQ S t'/Ti are constant of interaction of pair vise electron-scatterer U es (\q — Qj\ ja') and scatterer-scatterer 
U ss (\Qi — Qj \ /cr") potentials respectively, Qj are the scatterers positions (j = 1, ..N), a is the characteristic length 
of potential. 

Ensemble averaging on the configuration of interacting electron and classical heavy scatterers in quantum Kubo 
formulas has been performed according to the probability distribution proportional to \^\ , while the dynamic evolution 
has been realized according to recurrent relations (0) 
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The potential barrier of real type for electron -scatterer and scatterer-scatterer interaction has been taken in the 
Gaussian form: 

U(\q-q'\/a) = cxp(-\ q - q i\ 2 /a 2 ) (7) 
with equal to each other the all constants of interaction ( Q s = £g s > and a' = a"). 

V. NUMERICAL RESULTS 

Anderson localization of electrons have been investigated for one, two- and three dimensional disordered systems of 
scatterers. For ID case the results of our simulation and their brief description are presented below. In all calculations 
the dimesionless density of classical scatterer is approximately equal to unity (na ss 1), while ka ~ 1 and k~ 1 n ps 1. 

A. Electron conductivity of disordered systems of scatterers 

For the two small ratios of temperature to the height of a potential barrier (kT/Vo) equal to 0.04 and 0.0025 the 
Fig. |l| and Fig. || present ID results for diagonal elements of the tensor of the electrical conductivity, which are the 
momentum- momentum time correlation functions. Curves 1 on both figures relate to calculations taking into account 
only one term of iteration series (|5|), while curve 2 present results allowing for all terms of iteration series (^|). So 
curves 1 on both figures have been obtained by using only the classical trajectories without momentum jumps. The 
curve 1 on the Fig. [I] have practically the traditional fast decay. 

Curves 2 on both figures present momentum- momentum time correlation function obtained for dynamic trajectories 
with momentum jumps. The difference between these two correlation functions is larger at lower temperatures. 
Analyzing Fig. |^ one can conclude that curves 1 and 2 show undamped oscillations. Appearing of the undamped 
oscillations at very low temperature results from shortcomings of our model, in which the classical slow heavy scatterers 
are treated as non- moving particles on the time scale of the order of characteristic dynamic time of electrons. 

However as is known from the literature the undamped oscillations results in unphysical behaviour of the Fourier 
transform of the momentum- momentum time correlation function. The next Fig. [][]].[] and ^| show the significant 
difference of the Fourier transforms for correlation function with damped and undamped oscillations. Fig. [|, |i] and 
|5|, ^| demonstrate the real and imaginary parts of its Fourier transforms versus the dimensionless frequency Tluj/Vq. 
The real part of Fourier transform characterizes the Ohmic absorption of electromagnetic energy and has the physical 
meaning of electron conductivity, while the imaginary part presents (e — l)w, where e is permittivity of the system. 
The curve 2 (quantum trajectories) on Fig. || is higher than the curve 1 (classical trajectories) but both curves are 
going to zero at small frequency pointing out that the static conductivity at zero frequency is equal to zero or is 
very small. Note that points of curve 2 are going to zero faster than the same value obtained in approximation of 
classical trajectories (curve 1). This behaviour of conductivity in the vicinity of zero frequency is the characteristic 
manifestation of the electron Anderson localization. 

At lower temperatures the undamped oscillations with the proper phase shift of the momentum- momentum time 
correlation function may result in appearance of the negative values of real part of Fourier transform on Fig. ^. To 
overcome this shortcomings of our model at very low temperature it is necessary to take into account the slow motion 
of heavy particles, which should destroy the coherence oscillation of the light electrons trapped by heavy particles 
and cancel these negative values. 



B. Position and momentum dispersions 

The Fig. ^ presents the ID results for quantum position dispersion versus the dimesionless time tVo/H at the 
mentioned above two ratios of temperature to the height of a potential barrier (kT/Vo). Curves 1 and 2 correspond 
to kT/Vo equal to 0.04 and 0.0025 respectively. 

In one dimensional case quantum electrons is known Q should be localized at any temperatures but the length of 
localization is function of the particle energy. Analyzing data, submitted on the Fig. [?], it is possible to note that 
at lower temperature the curve 2 is very flat and the localization length (A) can be estimated as the squared root 
of the characteristic value of the position dispersion (A < 0.7). At higher temperature the behaviour of curve 1 is 
more complicated. At the initial stage (Wq/U < 100) curve 1 has a very flat part, which can be also considered as a 
characteristic manifestation of localization. 



G 



However when the time is more than 100 the curve 1 shows behaviour, which is typical for particle diffusion. This 
characteristic change in behaviour may be connected with momentum jumps of quantum trajectories. The virtual 
energy of the trajectories due to momentum jumps may be large. So the exponentially rare fast trajectories with very 
large value of localization length can give the exponentially large contribution to position dispersion as the difference 
in positions of fast trajectories may be exponentially large at large enough time. The similar problem connected with 
the exponentially large contribution of the exponentially rare events arises at the consideration of the classical wave 
propagation in random media. It is known Q for one dimensional case that the dispersion of wave intensity is not self- 
averaging value as the exponentially rare configurations of scatterers can give the exponentially large contribution of 
intensity at large enough distances from the wave source. Assuming that the initial flat part of curve 2 indicates the 
particle localization we can estimate the localization length as equal to 1.2. 

The Fig. |8| presents momentum dispersion of quantum particle in disordered systems of scatterers. Curves 1 and 
2 are stabilized and it is interesting to note that the lower temperature curve 2 are stabilized at larger value than 
curve 1. This may be connected with the uncertainty momentum- position principle and related fast oscillation of the 
localized quantum particles. 

C. Energy distribution function 

Quantum average energy of the system has been calculated according to the easily derived formula: 

If 1 

H{t) = 2v / dp 1 dq 1 dp 2 dq2-{H(p l7 q 1 ) + H(p2,q 2 )}x 

(zTrn) J z 



W (p 1 ,qi;p 2 ,q2;t;ihf3) 



Our calculations have shown that this function versus time variable is practically constant. To analyze the energy 
fluctuations during quantum dynamic evolution of the system the more interesting energy distribution function has 
been also calculated. The simple semi classical approximation of this function is defined by the following expression: 



p(E,0) = 

dt J dp 1 dq x dp2dq 2 \ i {8{E-H {pi,qi)) + S (E - H {p 2 , q 2 ))} x 



W (p 1 ,q 1 ;p2,q2;t;ihf3) 

The function In (p (E, (3)) versus E/Vq has been presented by curves 1 and 2 on the Fig. ^ and [l(] for higher and 
lower temperatures respectively. Curves 1 present results obtained in approximation of classical trajectories (only first 
term in iteration series), while curves 2 show results allowing for all terms of iteration series. At E/Vq < 4 curves 1 
and 2 practically coincide both at higher and lower temperatures. However at the E/Vq > 4 the distinction in curves 
1 and 2 is large due to momentum jumps resulting in appearance of trajectories with large virtual energy. So the 
asymptotic behaviour of the curves 1 and 2 is quite different. These asymptotics in logarithmic scale may be fitted by 
the strait lines with different slope. The existence of the long exponentially decreasing tail of the energy distribution 
function supports our explanation of the peculiarity in the behaviour of position dispersion of quantum particles. For 
detailed investigations the additional consideration is needed to analyze this effect. 

Lowering temperatures results also in spliting the main peak in energy distribution function into two sharper peaks. 
As it follows from Fig. |^ and [l^ the energy gap at the edge of the energy spectrum a bit increases with the lowering 
temperature. 

VI. CONCLUSION 

In the Wigner formulation of quantum statistical mechanics for canonical and micro canonical ensembles we have 
presented a new computational technique allowing quantum dynamics simulations for systems including subsystems of 
quantum interacting particles and subsystems of classical heavy scatterers as well as the system of quantum particles 
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in external potential field. The developed approach for quantum dynamics includes a sophisticated combination of 
well known molecular dynamics method and Monte Carlo technique. 

For electrons in disordered systems of scatterers the numerical results have been obtained for series of the average 
values of the quantum operators including position and momentum dispersions, average energy energy distribution 
function as well as for the frequency dependencies of tensor of electron conductivity and permittivity according to 
quantum Kubo formula. Zero or very small value of static conductivity have been considered as the manifestation of 
Anderson localization of electrons in ID case. Independent evidence of Anderson localization comes from the behaviour 
of the calculated time dependence of position dispersion. Nevertheless for localized electrons the energy distribution 
function obtained has the long exponentially decaying tail, which is the reason of exponentially rare appearance of 
large values of quantum particle virtual energy that strongly affects the behaviour of the position dispersion. 
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FIG. 1. Momentum-momentum time correlation function versus tVo/h for higher temperature (kT/Vo = 0.04): 1- approxi- 
mation of the classical trajectories; 2- quantum trajectories 

FIG. 2. Momentum-momentum time correlation function versus tVo/h for lower temperature (kT/Vo = 0.0025): 1- approx- 
imation of the classical trajectories; 2- quantum trajectories 

FIG. 3. The real part of Fourie transform of the momentum- momentum time correlation function versus hui/Vo for higher 
temperature (kT/Vo = 0.04): 1- approximation of the classical trajectories; 2- quantum trajectories 

FIG. 4. The real part of Fourie transform the momentum-momentum time correlation function versus Tuo/Vo for lower 
temperature (fcT/Vb=0.0025): 1- approximation of the classical trajectories; 2- quantum trajectories 

FIG. 5. The imaginary part of Fourie transform of the momentum-momentum time correlation function versus fkj/Vo for 
higher temperature (kT/Vo = 0.04): 1- approximation of the classical trajectories; 2- quantum trajectories 

FIG. 6. The imaginary part of Fourie transform the momentum-momentum time correlation function versus Hlu/Vo for 
lower temperature (fcT/Vb=0.0025): 1- approximation of the classical trajectories; 2- quantum trajectories 



FIG. 7. Position dispersion: 1- higher temperature (kTuj/Vo = 0.04); 2- lower temperature (kT/Vo = 0.0025) 
FIG. 8. Momentum dispersion: 1- higher temperature (kT/Vo = 0.04); 2- lower temperature (kT/Vo = 0.0025) 

FIG. 9. Energy distribution in logarithmic scale versus Eq/Vq for higher temperature (kT/Vo = 0.04): 1- approximation of 
the classical trajectories; 2- quantum trajectories 

FIG. 10. Energy distribution in logarithmic scale versus E /Vo for lower temperature (kT/Vo = 0.0025): 1- approximation 
of the classical trajectories; 2- quantum trajectories 
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